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December 13, 1888. 

Professor G. G. STOKES, D.C.L., President, in the Chair. 

The Presents received were laid on the table, and thanks ordered 
for them. 

The following Papers were read :— 

I. 44 Spectrum Analysis of Cadmium.” By A. Grunwald, 
Professor of Mathematics in the Imp. Roy. German Poly¬ 
technic University at Prague. Communicated by Professor 
Liyeing, F.R.S. Received November 26, 1888. 

[Publication deferred.] 


II. 44 On the Bending and Vibration of thin elastic Shells, espe¬ 
cially of Cylindrical Form.” By Lord Rayleigh, M.A., 
D.C.L., Sec. R.S. Received December 1, 1888. 

In a former publication* “ On the Infinitesimal Bending of Sur¬ 
faces of Revolution,” I have applied the theory of bending to explain 
the deformation and vibration of thin elastic shells, which are sym¬ 
metrical about an axis, and have worked out in detail the case where 
the shell is a portion of a sphere. The validity of this application 
depends entirely upon the principle that when the shell is thin 
enough and is vibrating in one of the graver possible modes, the 
middle surface behaves as if it were in extensible. “ When a thin 
sheet of matter is subjected to stress, the force which it opposes to 
extension is great in comparison with that which it opposes to bend¬ 
ing. Under ordinary circumstances, the deformation takes place 
approximately as if the sheet were inextensible as a whole, a condi¬ 
tion which, in a remarkable degree, facilitates calculation, though 
(it need scarcely be said) even bending implies an extension of all 
but the central layers.” If we fix our attention upon one of the 
terms involving sines or cosines of multiples of the longitude, into 
which, according to Fourier’s theorem, the whole deformation may 
be resolved, the condition of inextensibility is almost enough to 
define the type. If there are two edges, e.y., parallel to circles of lati- 

# ‘London Math. Soc. Proc./ vol. 13, p.-4, November, 1881. 
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tnde, the solution contains two arbitrary constants; but if a pole be 
included, as when the shell is in the form of a hemisphere, one of the 
constants vanishes, and the type of deformation is wholly determined, 
without regard to any other mechanical condition, to be satisfied at 
the edge or elsewhere. It will be convenient to restate, analytically, 
the type of deformation arrived at [equation (5)]. If the point 
upon the middle surface, whose coordinates were originally a, 9 , 0, 
moves to a + Sr, 0 + 80, 0 + £0, the solution is 

£0 = A tan s ijr 9 cos s(p q 

SO — — A sin 9 tan 5 i 9 sin $0 >.(1), 

Sr — A a (s-fcos 0) tan^ 9 sin s0 J 

9 being the colatitude measured from the pole through which the 
shell is complete. Any integral value higher than unity is admis¬ 
sible for s. The value 0 and 1 correspond to displacements not 
involving strain. 

In a recent paper* Mr. Love dissents from the general principle 
involved in the theory above briefly sketched, and rejects the special 
solutions founded upon it as inapplicable to the vibration of thin 
shells. The argument upon which I proceeded in my former paper, 
and which still seems to me valid, may be put thus : It is a general 
mechanical principlef that, if given displacements (not sufficient by 
themselves to determine the configuration) be produced in a system 
originally in equilibrium by forces of corresponding types, the result¬ 
ing deformation is determined by the condition that the potential 
energy of deformation shall be as small as possible. Apply this to an 
elastic shell, the given displacements being such as not of themselves 
to involve a stretching of the middle surface.J The resulting 
deformation will, in general, include both stretching and bending, 
and any expression for the energy will contain corresponding terms 
proportional to the first and third powers respectively of the thick¬ 
ness. This energy is to be as small as possible. Hence, when the 
thickness is diminished without limit, the actual displacement will be 
one of pure bending, if such there be, consistent with the given con¬ 
ditions. Otherwise the energy would be of the first order (in 
thickness) instead of, as it might be, of the third order, in violation 
of the principle. 

It will be seen that this argument takes no account of special 
conditions to be satisfied at the edge of the shell. This is the point 
at which Mr. Love concentrates his objections. He considers that 

# “ On the small free Yibrations and Deformation of a thin elastic Shell,” ‘ Phil. 
Trans.,’ A, 1888. 

f i Phil. Mag.,’ March, 1875 ; * Theory of Sound,’ § 74. 

J There are cases where no displacement (involving strain at all) is possible 
without stretching of the middle surface, e.g., that of the complete sphere. 
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the general condition necessary to be satisfied at a free edge is in 
fact violated by such a deformation as (1). But the condition in 
question* contains terms proportional to the first and to the third 
powers respectively of the thickness, the coefficients of the former 
involving as factors the extensions and shear of the middle surface. 
It appears to me that when the thickness is diminished without limit, 
the fulfilment of the boundary condition requires only that the middle 
surface be unstretched, precisely the requirement satisfied by solutions 
such as (1). 

Of course, so long as the thickness is finite, the forces in operation 
will entail some stretching of the middle surface, and the amount of 
this stretching will depend on circumstances. A good example is 
afforded by a circular cylinder with plane edges perpendicular to the 
axis. Let normal forces locally applied at the extremities of one 
diameter of the central section cause a given shortening of that 
diameter. That the potential energy may be a minimum, the de¬ 
formation must assume more and more the character of mere bending 
as the thickness is reduced. The only kind of bending that can occur 
in this case is the purely cylindrical one in which every normal 
section is similarly deformed, and then the potential energy is propor¬ 
tional to the total length of the cylinder. We see, therefore, that if 
the cylinder be very long, the energy of bending corresponding to the 
given local contraction of the central diameter may become very 
great, and a heavy strain is thrown upon the principle that the 
deformation of minimum energy is one of pure bending. 

If the small thickness of the shell be regarded as given, a point 
will at last be attained when the energy can be made least by a 
sensible local stretching of the middle surface such as will dispense 
with the uniform bending otherwise necessary over so great a length. 
But even in this extreme case it seems correct to'say that, when the 
thickness is sufficiently reduced, the deformation tends to become one 
of pure bending. 

At first sight it may appear strange that of two terms in an ex¬ 
pression of the potential energy, the one proportional to the cube of 
the thickness is to be retained, while that proportional to the first 
power may be omitted. The fact, however, is that the large potential 
energy which would accompany any stretching of the middle surface 
is the very reason why such stretching will not occur. The compara¬ 
tive largeness of the coefficient (proportional to the first power of the 
thickness) is more than neutralised by the smallness of the stretching 
itself, to the square of which the energy is proportional. 

In general, if Yq be the coordinate measuring the violation of the 
tie which is supposed to be more and more insisted upon by increasing 


See liis equation (33). 
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stiffness, and if the other coordinates be suitably chosen, the potential 
energy of the system may be expressed 

V = ic 1 Y' 1 3 +|c 2 Y'2 3 + k c sfi+ . 

This follows from the general theorem that Y and T may 
always be reduced to sums of squares simply, if we suppose that 
T — 

The equations of equilibrium under the action of external forces 
Sfq, ^ 2 , .. are thus 

*1 = C lfv *2 = C 2^2).; 

hence if the forces are regarded as given, the effect of increasing c l 
without limit is not merely to annul y^, but also the term in V which 
depends upon it. 

An example might be taken from the case of a rod clamped at one 
end A, and deflected by a lateral force, whose stiffness from the end 
A up to a neighbouring place B, is conceived to increase indefinitely. 
In the limit we may regard the rod as clamped at B, and neglect the 
energy of the part AB, in spite of, or rather in consequence of, its 
infinite stiffness. 

If it be admitted that the deformations to be considered are pure 
bendings, the next step is the calculation of the potential energy 
corresponding thereto. In my former paper, the only case for which 
this part of the problem was attempted was that of the sphere. 
A fter bending, “ the principal curvatures differ from the original 
curvature of the sphere in opposite directions, and to an equal amount,* 
and the potential energy of bending corresponding to any element of 
the surface is proportional to the square of this excess or defect of 
curvature, without regard to the direction of the principal planes.’* 
Though he agrees with my conclusions, Mr. Love appears to regard 
the argument as insufficient. But clearly in the case of a given 
spherical shell, there are no other elements upon which the energy of 
bending could depend. “Thus the energy corresponding to the 
element of surface u 2 sind d6 d<p may be denoted by 

a 3 H(£/> -1 ) 3 sin0 d6 dcj) . . (2), 

where H depends upon the material and upon the thickness.” 

By the nature of the case H is proportional to the elastic constants 
and to the cube of the thickness, from which it follows by the method 
of dimensions that it is independent of a, the radius of the sphere. 

* This is in virtue of Gauss’s theorem that the product of the principal curva¬ 
tures is unaffected by bending. 
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I did not, at the time, attempt the further determination of H, not 
needing it for my immediate purpose. Mr. Love has shown that 

H = | -nh?> .(3), 


where 2 h represents the thickness, and n is the constant of rigidity. 
Why n alone should occur, to the exclusion of the constant of com¬ 
pressibility, will presently appear more clearly. 

The application of (2) to the displacements expressed in (1) gave 
[equation (18)] 

Y = 2 (s 3 — s) A s 2 H sin -8 # tan 2s ^0d0 .(4), 

Jo 


0 being the colatitude of the (circular) edge. In the case of the 
hemisphere of uniform thickness 

Y = (s 3 -s) (2s 2 - 1) A s 2 . (5). 

The calculation of the pitch of free vibration then presented no 
difficulty. If a denote the superficial density, and cos pt represent the 
type of vibration, corresponding to s == 2, p s to s = 3, and so on, it 
appeared that 


v 2 = -^x 5-2400, 


Ts = 


L? x 14726, 

n.*cr 


Pi = 


v/H 


x 28-462 


so that 

Ps/Pa ~ 2-8102, pjp s = 5-4316, 

determining the intervals between the graver notes. 

If the form of the shell be other than spherical, the middle surface 
is no longer symmetrical with respect to the normal at any point, and 
the expression of the potential energy is more complicated. The 
question is now not merely one of the curvature of the deformed 
surface; account must also be taken of the correspondence of normal 
sections before and after deformation.* A complete investigation 
has been given by Love; but the treatment of the question now to be 
explained, even if less rigorous, may help to throw light upon this 
somewhat difficult subject. 

In the actual deformation of a material sheet of finite extent there 
will usually be at any point not merely a displacement of the point 
itself, but a rotation of the neighbouring parts of the sheet, such as a 

# An extreme case may serve as an illustration. Suppose that the bending is 
such that the principal planes retain their positions relatively to the material sur¬ 
face, but that the principal curvatures are exchanged. The nature of the curvature 
at the point in question is the same after deformation as before, and by a rotation 
through 90° round the normal the surfaces may be made to lit; nevertheless the 
energy of bending is finite. 
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rigid body may undergo. All this contributes nothing to the energy. 
In order to take the question in its simplest form, let us refer the 
original surface to the normal and principal tangents at the point in 
question as axes of coordinates, and let us suppose that after deforma¬ 
tion, the lines in the sheet originally coincident with the principal 
tangents are brought back (if necessary) to occupy the same positions 
as at first. The possibility of this will be apparent when it is re¬ 
membered that in virtue of the in extensibility of the sheet, the angles 
of intersection of all lines traced upon it remain unaltered. The 
equation of the original surface in the neighbourhood of the point 
being 



that of the deformed surface may be written 


z 


-x f. 

— 2 


1 TT~ + —t * 7X V ( 
iPl + fy 1 P%+°P*i J 


4- 2 T-xy 


\ 


( 7 ). 


In strictness (p l + °Pi)~ l ,(P 2 4- are the curvatures of thesections 

made by the planes x = 0, y = 0; but since principal curvatures are 
a maximum or a minimum, they represent with sufficient accuracy 
the new principal curvatures, although these are to be found in 
slightly different planes. The condition of inextensibility shows that 
points which have the same x and y in (6) and (7) are corresponding 
points, and by Gauss’s theorem it is further necessary that 


?& + *>:= 0 . . . ( 8 ). 

Pi Pst ... 

It thus appears that the energy of bending will depend upon two 
quantities, one giving the alterations of principal curvature, and the 
other t depending upon the shift (in the material) of the principal 
planes. 

In calculating the energy we may regard it as due to the stretchings 
and contractions under tangential forces of the various infinitely thin 
laminse into which the shell may be divided. The middle lamina 
being unstretched, makes no contribution. Of the other lamines, the 
stretching is in proportion to the distance from the middle surface, 
and the energy of stretching is therefore as the square of this 
distance. When the integration over the whole thickness of the shell 
is carried out, the result is accordingly proportional to the cube of 
the thickness. 

The next step is to estimate more precisely the energy correspond¬ 
ing to a small element of area of a lamina. The general equations in 
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three dimensions, as given in Thomson and Tait’s * Natural Philosophy,’ 
§ 694, are 



na = S, 

nb = T, nc =U .. , 

.( 9 ) 



Me = P — a (Q + R)'| 

M/ = Q —a (R+P) V. 

M^ = R-<7(P + Q)J 

....(10), 

where 


1 1 
§ 

II 

<5. 

...(H).* 


The energy w, corresponding to the unit of volume, is given by 
2w = (m + n ) (h 2 +/ 3 -f- 

4-2 (m — n) (fg+ge f e/)-f n(a 2 + 6 3 4c 2 ).. (12). 

In the application to a lamina, supposed parallel to xy, we are to 
take It — 0, S = 0, T = 0; so that 

e+/ A 7 A 

g — —G --, a = 0 , & = 0 . 

1 —a 

Thus in terms of the elongations e , f parallel to x , y, and of the 
shear c, we get 


« = «{rf»+/»+^ (*+/)•+iri*}.(13). 

We have now to express the elongations of the various laminae of a 
shell when bent, and we will begin with the case where t = 0, that is, 
when the principal planes of curvature remain unchanged. It is 
evident that in this case the shear c vanishes, and we have to deal 
only with the elongations e and f parallel to the axes. In the section 
by the plane of zx, let s , s' denote corresponding infinitely small 
arcs of the middle surface and of a lamina distant h from it. If ^ 
be the angle between the terminal normals, s = p^, s' = ( p 1 4 K) yfr, 
s' — s == h^jr. In the bending, which leaves s unchanged, 

S$' = hhyjr = hsh (1 jp{). 

Hence 

e = Ss/s' = hS^l/p^, 

and in like manner / & (l/p 2 )* Thus for the energy U per unit of 

area we have 

* M is Young’s modulus, cr is Poisson’s ratio, n is the constant of rigidity, 
and {m — \n) that of cubic compressibility. In terms of Lame’s constants (a., /x), 
m = A + p, n = fi. 
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dJJ = nip 1 dli 


IA Pi/ \ /V + p x p % J J 


and on integration oyer the whole thickness of the shell (2 h) * 
2 nh s 


U = ■ 


, {( i 7;T + (^T + s+-:(^ +! 7;T}- <M) ' 


This conclusion may he applied afc once, so as to give the result 
applicable to a spherical shell; for, since the original principal planes 
are arbitrary, they can be taken so as to coincide with the principal 
planes after bending. Thus r = 0; and by Gauss’s theorem 


so that 


1 1 

a—■= o, 

Pi Pi 



(15), 


where Bp~ l denotes the change of principal curvature. Since e = —/, 
g — 0, the various laminse are simply sheared, and that in proportion 
to their distance from the middle surface. The energy is thus a 
function of the constant of rigidity only. 

The result (14) is applicable directly to the plane plate; but this 
case is peculiar in that, on account of the infinitude of p v Pn (8) is 
satisfied without any relation between fyq and dp 2 . Thus for a plane 
plate 


tt _ f 11 m—n f 1 1 VI 

3 L pi 2 pJ* m+n\p l Pi) 1 


(16), 


where p x l , p%~ 1 , are the two independent principal curvatures after 
bending. 

We have thus far considered t to vanish ; and it remains to inves¬ 
tigate the effect of the deformations expressed by 


$z = rxy = lT(f 3 — r?) ...... ...(17), 

where y relate to new axes inclined at 45° to those of x , y. The 
curvatures defined by (17) are in the planes of J, y, equal in numerical 
value and opposite in sign. The elongations in these directions for 

* It is here assumed that m and n are independent of h> that is, that the material 
is homogeneous. If we discard this restriction, we may form the conception of a 
shell of given thickness, whose middle surface is physically inextensible, while yet 
the resistance to bending is moderate. In this way we may realise the types of 
deformation discussed in the present paper, without supposing the thickness to he 
infinitely small; and the independence of such types upon conditions to he satisfied 
at a free edge is perhaps rendered more apparent. 
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any lamina within the thickness of the shell are hr , — hr, and the 
corresponding energy (as in the case of the sphere just considered) 
takes the form 


IT = 


4»7i s T® 

3 


.(18). 


This energy is to be added* to that already found in (14) ; and we 
get finally 


U = 


2 f B 
3 


{ (°7, + 8 sT + 2t ’} • ■ (m 


is the complete expression of the energy, when the deformation is 
Mich that the middle surface is unextended. We may interpret r by 
nans of the angle x f through which the principal planes are shifted ; 
>hus 



It will now be in our power to treat more completely a problem of 
great interest, viz., the deformation and vibration of a cylindrical 
shell. In my former paper I investigated the types of bending, but 
without a calculation of the corresponding energy. The results were 
as follows.f If the cylinder be referred to columnar coordinates 
z, r, 0, so that the displacements of a point whose equilibrium co¬ 
ordinates are z, a, 0 are denoted by hz, hr, ah(p, the equations express¬ 
ing inextensibility take the form 


** 1 « 

II 

o 

_ dh(j) 

Sr+a !$ = 0 ’ 

2 dc(J> 

__ +a2 __ = 0 ... 

.. (21), 

from which we may 

deduce 




dPty _ 
dz* 


(22). 


By (22), if <50 cc cos s0, we may take 


ah<p = (A*a+B* 0 ) cos s<fi ................ (23), 

and then, by (21) 

hr = s (A,a4-B s z) sin 50 ......... .... (24), 

hz = — s^B^a sin s0 ................. (25). 


# There are clearly no terras involving the products of r with the changes of 
principal curvature ^(Ps -2 ) 5 f° r a change in the sign of r can have no 

influence upon the energy of the deformation defined by (7). 

f The method of investigation is similar to that employed by Tellet in his 
memoir ( <£ On the Properties of Inextensible Surfaces,” * Irish Acad. Trans.,’ voh 22, 
1855, p. 179), to which reference should have been made. 
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If the cylinder be complete, s is integral; A s and are independent 
constants, either of which may vanish. In the latter case the dis¬ 
placement is in two dimensions only.* It is unnecessary to stop to 
consider the demonstrations of (21), inasmuch as these equations 
will present themselves independently in the course of the investiga¬ 
tions which follows. 

It will be convenient to replace 8z, £r, n£0 by single letters, which, 
however, it is difficult to choose so as not to violate some of the usual 
conventions. In conformity with Mr. Love’s general notation, I will 
write 

hz — u, ahcf) — v, — iv .. ( 26 ). 

The problem before us is the expression of the changes of principal 
curvature and shifts of principal planes at any point P (z, 0) of the 
cylinder in terms of the displacements u, v , w. As in (6), take as 
fixed co-ordinate axes the principal tangents and normal to the 
undisturbed cylinder at the point P, the axis of x being parallel to 
that of the cylinder, that of y tangential to the circular section, and 
that of f normal, measured inwards. If, as it will be convenient to 
do, we measure z and 0 from the point P, we may express the undis¬ 
turbed coordinates of a material point Q in the neighbourhood of P, 
by 

x = z, y — f . (27). 

During the displacement the coordinates of Q will receive the 
increments 

u, } w sin 0 -f v cos 0, —w cos 0 + v sin 0 ; 
so that after displacement 

x = z + u, y — a0-fu’0-J-?;(l— ^0 2 ), 

f — w(l — 

or if u , v, w be expanded in powers of the small quantities z, 0, 


du du , , 

w — —2 ~h "TT0 4- .. 

ciZq ct(p © 

.(28). 

, . dv dv . 

y = a0+« <( 0+« o + 5 - a +^+. ...... 

.(29). 


* See 4 Theory of Sound,’ § 233. 
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„ i dm dw 

£ = j;a<t?-w 0 --^ z-^f+v o0 


d<l >o 




dz 0 drf>Q Z< ^ 2 d(f) 0 2 


, dv , , dv , 9 
+ “*0+—0* 


dZn 


#0 


■ (30), 


v 0 , . . . being the values of u, v at the point P. 

These equations give the coordinates of the various points of the 
deformed sheet. We have now to suppose the sheet moved as a rigid 
bodj so as to restore the position (as far as the first power of small 
quantities is concerned) of points infinitely near P. A purely translatory 
motion by which the displaced P is brought back to its original posi¬ 
tion will be expressed by the simple omission in (28), (29), (30) of 
the terms u 0 , v 0 , iv 0 respectively, which are independent of z, 0. The 
effect of an arbitrary rotation is represented by the additions to 
x, y , f respectively of yO s — g*<9 3 , <£6 X — x6 s , x0 2 — y0 r ; where for the 
present purpose 0 1? 0 3 are small quantities of the order of the 

deformation, the square of which is to be neglected throughout. If 
we make these additions to (28), &c., substituting for x, y, g* in the 
terms containing 0 their approximate values, we find so far as the 
first powers of 0 


dv dv . ' 

y = a 9 +w 0 <p+j-z+^-zO s , 


div dw , , . 

£■= — 9 i + v o 0+20 3 -a00 1 . 


#0 


How, since the sheet is assumed to be inextensible, it must be pos¬ 
sible so to determine ®\-> % ^3 that to this order x — z,y~ a0, g' = 0. 


Hence 

du 

dz 0 

o, 


du 

■xr+ a0 3 : 

d4> 0 

= 0 , 

du 


0 , 


dv 

= 0, 

dz 0 

-*3 = 


w o + TT 
d<P 0 

dw 


0 , 

dw 

— v 0 +ae 1 

= 0. 

dz 0 

+ 

II 

d*fi 0 


The conditions of inextensibility are thus (if we drop the suffices 
as no longer required) 

VOL. xlv. 


1 
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w+ ^ =0 > 


du dv 

_ + a _ = 0 ... (3i), 


which agree with (21), 

Returning to (28), &c., as modified by the addition of the trans- 
latory and rotatory terms, we get 


a? = 0 -f terms of 2nd order in z, 0, 

y = O0+ 




r n m 




#<f 


dz^clcpQ 


00 


2 d0 o 3 cfe 0 ^ rf0 o 1 


or since by (31) d^wjdz 9, = 0, and dvjd(j> — —w> 
£ — ia(j) 2 —iw o 0 3 


d%w . , dhet , dv , 

'<k Q d<f^ * 11 %^ + ir 0 ^' 


The equation of the deformed surface after transference is thus 

„ _ r i *> i* i , ,ji i i 

'■ “ i a <fe 0 a tfz o d0 o J + ^ 12a 2a 2 W ° 2a 2 # 0 2 j ‘ * * * ^ 32) ' 

Comparing with (7) we see that 

S—=0, 

Pi 

so that by (19) 

tr = — ( _£L 1 f „+£» V+ C 

3a 2 lm+» a® V + j + V 



1 / dhv \ 

_ _1 /dv 

d^w \ 

5—- 
P* ~ 


a 

dzd<p) 


du__d%\2-l 
.<& dzd<t>) j ' ’ * l 


This is the potential energy of bending reckoned per unit of area. 
It can if desired be expressed by (31) entirely in terms of v* 

We will now apply (24) to calculate the whole potential energy of 
a complete cylinder, bounded by plane edges z = +Z, and of thick* 

* From Mr. Love’s general equations (12), (13), (18) a concordant result maybe 
obtained by introduction of tbe special conditions— 

h ~ 0, Ji» — lja, 1 /pi — 0, l/p 2 = lja } 

limiting tbe problem to tbe case of tbe cylinder, and of those 

<r 2 t= ff 2 =. m = 0, 

which express the inextensibility of tbe middle surface. 
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ness which, if variable at all, is a function of z only. Since u, v , w 
are periodic when 0 increases by 27r, their most general expression in 
accordance with (31) is [compare (23), &c.] 


v = 2[(A 5 a + B,, 2 ) cos 50—(A/a + B/z) sin 50]. (35), 

w = 2 [5 (A 5 a-f B^z) sins0 + 5 (A/a+B/ 2 ) cos 50],... (36), 
u = 2[ —s^B^a sin 50—cos 50],.. (3?), 


in which the summation extends to all integral values of s from 0 
to co. But the displacements corresponding to s == 0, s = 1 are 
such as a rigid body might undergo, and involve no absorption of 
energy. When the values of u , v, w are substituted in (34) all the 
terms containing products of sines or cosines with different values of 
s vanish in the integration with respect to 0, as do also those which 
contain cos 50 sin 50. Accordingly 


|- 2 tt 

I Uad0 = 
Jo 


47 mh? 

3 a 


r m 1 

L m+n a 3 


( s 3 -*) 3 


{ (A,a+B,«? + (A;a+B/*)*}+2 (*-!)»(B/+B/’)J .. (38). 


Thus far we might consider k to be a function of z ; but we will no w 
treat it as a constant. In the integration with respect to z the odd 
powers of z will disappear, and we get as the energy of the whole 
cylinder of radius a , length 2 1, and thickness 2 h, 



U a d(j)dz 


87 mliH 
3a 2 



A/ + A/ 3 


+ oj (B/+B/2)} + B» s + B;»].(39), 

in which s = 2, 3, 4, ... . 

The expression (39) for the potential energy suffices for the solu¬ 
tion of statical problems. As an example we will suppose that the 
cylinder is compressed along a diameter by equal forces F, applied at 
the points z = z v 0 = 0, 0 = «r, although it is true that so highly 
localised a force hardly comes within the scope of the investigation in 
consequence of the stretchings of the middle surface, which will 

1 2 
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dccur in the immediate neighbourhood of the points of applica¬ 
tion.^ 

The work done upon the cylinder by the forces F during’ the 
hypothetical displacement indicated by dA s , &c., will be by (86) 

—F2 s (cidA s ' (1 + cos stt), 

so that the equations of equilibrium are 


0 

II 

^ 1^3 

© 

II 

•§|g 

—t—7 — — (1 + C03 S7r) saF. 

dA $ ' v ' 

= — (1 -f COS S7r) SZ] F 

Cuxj s 


Thus for all values of s , 

A s = B = 0; 


and for odd values of s , 


a; = B/ = 0. 


But when s is even, 


ws 3 . ,_ 3sn 2 F 

m + n s 87rnhH(s% —l) 2 

f Z 2 1 __ Ssaz{F 

\m-\-n 8a 2 J 8 87 mh s l(s^ —l) 2 

and the displacement it? at any point (z, 0) is given by 

w = 2(A 2 'a~}-B 2 V) cos 20+ 4 (A/a+B/ 2 ) cos 40+ ... 


(40), 


( 41 ); 


(42), 


where A 2 , B 2 ', A 4 f , .... are determined by (40), (41). 

If the cylinder be moderately long in proportion to its diameter, 
the second term in the left hand member of (41) may be neglected, so 
that 

A/ 

8a 2 z l a 


In this case (42) may be written 


= ^1 + { 2 A % a cos 20 + 4A l a cos 40 +. . . .}... (43), 


# Whatever the curvature of the surface, an area upon it may he taken so small 
as to behave like a plane, and therefore bend, in violation of Gauss’s condition, when 
subjected to a force which is so nearly discontinuous that it varies sensibly within 
the area. 
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showing that, except as to magnitude and sign, the curve of deforma¬ 
tion is the same for all values of z 1 and a* 

If 3 = + the amplitudes are in the ratio 1 + 3 zfjl 2 ; and if, 
further, z 1 = l , i.e., if the force he applied at one of the ends of the 
cylinder, the amplitudes are as 2 : —1. The section where the 
deformation (as represented by to) is zero, is given by 3= 0, 
in which if z 1 ■= l, z = —JL 

When the condition as to the length of the cylinder is not imposed, 
the ratio B/ : A/ is dependent upon s, and therefore the curves of 
deformation vary with z 9 apart from mere magnitude and sign. If, 
however, we limit ourselves to the more important term 5 = 2, we 
have 

4m AJ f 4m l 2 1 B 2 ' 

—-— <j ————^ -j-1 >-> 

m-\-n a ..... L m-hn oar J z x 

and » = 2B 2 ' { ~(~ + +* } cos 20; 

so that w vanishes when 


P ^m + n 
a 2 Sa 2 4 m 


(44). 


This equation may be applied to find what is the length of the 
cylinder when the deformation just vanishes at one end if the force 
is applied at the other. If z l = —* z = Z, 

l_ __ / f 3 (m+n) 1 

a \/ l 8772, j 

For many materials a [equation (11)] is about or m = 2 n. In 
such cases the condition is 


l &a. 

It should not be overlooked that although iv may vanish, u remains 
finite. 

Reverting to (23), (24), (25) we see that, if the cylinder is open at 
both ends, there are two types of deformation possible for each value 
of s. If we suppose the cylinder to be closed at z = 0 by a flat disk 
attached to it round the circumference, the inextensibility of the disk 
imposes the conditions, w = dr = 0, v = ah(j) = 0, when 3 = O.f 
Hence A $ = 0, and the only deformation now possible is 

* That w is unaltered when z and z l are interchanged is an example of the general 
law of reciprocity. 

f s being greater than 1. 
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w — cr = stt s z sm 50 J 
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Another disk, attached where z has a finite value, would render the 
cylinder rigid. 

Instead of a plane disk let us next suppose that the cylinder is 
closed at z. = 0 by a hemisphere attached to it round the circum¬ 
ference. By (1) the three component displacements at the edge of the 
hemisphere (0 = r) are of the form 

v = a S(f> = a cos 50. 
u = ctcO = — a sin 50. 
w = 8r = sa sin 50. 

Equating these to the corresponding 1 values for the cylinder, as given 
hy (23), (24), (25), we get 

A $ 1, B,y ^ S J 

so that the deformation of the cylinder is now limited to the type 
v = (a+sz) cos 50 

w — s (a+sz) sin 50 )> .(46), 

u = —a sin 50 

in which we may, of course, introduce an arbitrary multiplier and an 
arbitrary addition to 0. If the convexity of the hemisphere be turned 
outwards, z is to be considered positive. 

In like manner any other convex additions at one end of the 
cylinder might be treated. There are apparently three conditions to 
be satisfied by only two constants, but one condition is really re¬ 
dundant, being already secured by the inextensibility of the edg’es 
provided for in the types of deformations determined separately for 
the two shells. Convex additions, closing both ends of the cylinder, 
render it rigid, in accordance with Jellet’s theorem that a closed 
oval shell cannot be bent. 

It is of importance to notice how a cylinder, or a portion of a 
cylinder, can not be bent. Take, for example, an elongated strip, 
bounded by two generating lines subtending at the axis a small 
angle. Equations (31) [giving dhvjdr* = 0] show that the strip 
cannot be bent in the plane containing the axis and the middle 
generating line.* The only bending symmetrical with respect to this 

* This is the principle upon which metal is corrugated. 
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plane is a purely cylindrical one which leaves the middle generating 
line straight. There are two ways in which we may conceive the 
strip altered so as to render it susceptible of the desired kind of bend¬ 
ing. The first is to take out the original cylindrical curvature, which 
reduces it to a plane strip. The second is to replace it by one in which 
the middle line is curved from the beginning, like the equator of a 
sphere or ellipsoid of revolution. In this case the total curvature 
being finite, the Gaussian condition can be satisfied by a change of 
meridianal curvature compensating the supposed change of equatorial 
curvature. It is easy to calculate the actual stiffness from (8) and 
(14), for here t = 0. We have 


U = 


2 nlfi 
3 



1 1 pf , 

p<f m + n 



(47), 


which expresses the work per unit of area corresponding to a given 
bending $Pi~ l along the equator. If p x = go, the cylindrical strip is 
infinitely stiff. If the curvature be spherical, /> 2 — p l5 and 


and if = oc, 



Whatever the equatorial curvature may be, the ratio of stiffnesses in 
the two cases is equal to m : m + w, or about 2 : 3, the spherically 
curved strip being the stiffer. 

The same principle applies to the explanation of Bourdon’s gauge. 
In this instrument there is a tube whose axis lies along an arc of a 
circle and whose section is elliptical, the longer axis of the ellipse 
being perpendicular to the general plane of the tube. If we now 
consider the curvature at points which lie upon the axial section, we 
learn from Gauss’s theorem that a diminished curvature along the 
axis will be accompanied by a nearer approach to a circular section, 
and reciprocally. Since a circular form has the largest area for a 
given perimeter, internal pressure tends to diminish the eccentricity 
of the elliptic section and with it the general curvature of the tube. 
Thus, if one end be fixed, a pointer connected with the free end may 
be made to indicate the internal pressure.^ 

* Dec. 19.—It appears, however, that the bending of a curved tube of elliptical 
section cannot be pure, since the parts of the walls which lie furthest from the 
circular axis are necessarily stretched. The difficulty thus arising may be obviated 
by replacing the two halves of the ellipse, which lie on either side of the major 
axis, by two symmetrical curves which meet on the major axis at a finite angle . 
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W e will now proceed with the calculation for the frequencies of 
vibration of the complete cylindrical shell of length 21. If the 
volume-density be p* we have as the expression of the kinetic energy 
by means of (35), (36), (37). 

T = \ . 2 hp . Jj* (tt 2 + ^ + w 3 ) adcfrdz 
= 2 r n‘pdfhld S {a 3 (l-fs 2 ) (A,+A/ 3 ) 


+ [^(i +s 2 ) + s - V ] (B,« + B/»)}.(50). 

From these expressions for Y and T in (39), (50) the types and fre¬ 
quencies of vibration can be at once deduced. The fact that the 
squares, and not the products, of A s , B s , are involved, shows that 
these quantities are really the principal coordinates of the vibrating 
system. If A s , or A,', vary as cos jo s t, we have 




mn h? (s 3 — s) 3 
i m + npa 4 s 3 +1 


(51). 


This is the equation for the frequencies of vibration in two dimen¬ 
sions, f For a given material, the frequency is proportional to the 
thickness and inversely as the square on the diameter of the 
cylinder. $ 

In like manner if B s , or B/, vary £ as cos we find 


3a 3 m + n 

, 2 4 mn k 3 (s 3 —s) 3 s 3 l 3 m 

3 m+n pa 4 s 3 + I 3a 3 

+ (s 4 + s 3 )Z 3 


(52). 


If the cylinder be at all long in proportion to its diameter, the 


According to the equations (in columnar co-ordinates) of my former paper, the 
conditions that dr, dz shall be independent of (p lead to— 


dr = Cr, 




= 0 , 


where C is an absolute constant. 

The case where the section is a rhombus (drjdz = ± tan a) may be mentioned. 

The difficulty referred to above arises when drjdz — oo. 

* This can scarcely be confused with the notation for the curvature in the pre= 
ceding parts of the investigation. 

f See ‘ Theory of Sound,’ § 233. 

X There is nothing in these laws special to the cylinder. In the case of similar 
shells of any form, vibrating by pure bending, the frequency will be as the thick¬ 
nesses and inversely as corresponding areas. If the similarity extend also to the 
thickness, then the frequency is inversely as the linear dimension, in accordance 
with the general law of Cauchy. 
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difference between jp s ' and p s becomes very small. Approximately in 
tbis case 


Palp* — 14 - 


oa 2 /m + n 1 \ # 


or if we take m = 2 w, 5 


2 , 

Pa Ip 2 = 1 + 2ol2 • 


In my former paper I gave the types of vibration for a circular 
cone, of which the cylinder may be regarded as a particular case. 


In terms of columnar coordinates (z, r, 0 ) we have 

£0 = cos 50 .. (53), 

Sr = 5 tan 7 (A^-f B*) sin 50 . (54), 

Sz = tan 2 7 s (A^d-B*)] sin 50 . (55)* 


7 being the semi-vertical angle of the cone. Bor the calcula¬ 
tion of the energy of bending it would be simpler to use polar co¬ 
ordinates ( r , d, 0 ), r being measured from the vertex instead of from 
the axis. 

If the cone be complete up to the vertex, we must suppose, in 
(53) &c., B 5 = 0 . And if we proceed to calculate the potential 
energy, we shall find it infinite, at least when the thickness is uni¬ 
form. Bor since A s is of no dimensions in length, the square of the 
change of curvature must be proportional to A s 2 z~%. When this is 
multiplied by zdz, and integrated, a logarithm is introduced, which 
assumes an infinite value when z = 0. The complete cone must 
therefore be regarded as infinitely stiff, just as the cylinder would be 
if one rim were held fast. 

If two similar cones (bounded by circular rims) are attached so 
that the common rim is a plane of symmetry, the bending may be 
such that the common rim remains plane. If the distance of this plane 
from the vertex be Zj, the condition to be satisfied in (53) &c., is 
that Sz = 0 where z = z 2 . Hence 

g 0 = A s 11 — ~ | cos 50 .(56), 

sin 50 .(57)? 

Sz — 5 tan 2 7 A s {zj—z} sin 50 . (58), 


Sr ■ 


? tan 7 A s (z } 









